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SUMMARY 


The  problem  of  pressure  fluctuations  at  a  rigid  wall  under  a  turbulent 
boundary  layer  has  attracted  much  attention  in  the  past  decade.  At  low 
Mach  numbers  the  theory  la  well  established  from  the  work  of  Kralohnan 
and  Lllley,  and  reasonable  agreement  la  obtained  with  the  experiments  of 
Wlllmarth,  Hodgson  and  others.  At  high  Mach  numbers,  measurements  exist 
due  to  the  work  of  Klstler  and  Chen  but  so  far  no  theory  Is  available, 
apart  from  that  due  to  Phillips,  whloh  Is  however  related  to  the  noise 
radiated  from  supersonic  turbulent  shear  flows. 

This  Report  reviews  the  theory  of  wall  pressure  fluctuations  in 
incompressible  flow,  and  shows  how  tho  character  of  the  pressure  fluc¬ 
tuations  changes  In  passing  from  the  flow  to  the  wall,  Attention  Is 
drawn  to  the  more  Important  interactions  giving  rise,  to  the  pressure 
fluotuatlona,  as  well  as  to  the  region  of  the  boundary  layer  mainly  res¬ 
ponsible  for  the  wall  pressure  fluctuations, 

The  work  la  extended  to  include  the  effeota  of  oompreaalblllty.  It  la 
found  that  an  analysis  similar  to  that  of  Phillips  la  appropriate,  although, 
unlike  the  latter  work,  this  new  treatment  la  not  reatrloted  to  the  oaae 
of  very  high  supersonic  Maoh  numbers,  The  analyala  makes  use  of  the  ratio 
aw/ur  ns  a  largo  parameter,  where  a„  Is  the  speed  of  sound  at  the  wall 
and  ur  la  the  shear  veloolty.  This  Is  certainly  true  for  a  very  wide 
range  of  Maoh  numbers  provided  that  the  wall  la  not  subjected  to  large 
rates  of  heat  transfer,  It  la  shown  that  the  wall  pressure  fluotuatlona 
are  now  the  result  of  fluctuations  In  both  the  vortlolty  snd  sound  modes, 

At  high  Maoh  numbers,  the  latter  contribution  is  In  the  form  or  eddy  Maoh 
waves,  as  suggested  by  Phillips.  On  making  certain  assumptions  regarding 
the  dominant  interactions,  estimates  of  the  magnitude  and  spectrum  of  the 
well  pressure  fluctuations  are  made  which  show  similar  trends  to  tho 
measurements  of  Klstler  and  Chon, 
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SOMMAIRE 


Lo  problems  des  fluotuations  de  presalon  a  une  parol  rlglde  sous  une 
couohe  limitrophe  turbulente  a  souleve  une  grande  attention  durant  oette 
derntlre  dloade,  A  dee  nombres  de  Maoh  peu  sieves,  la  theorie  eat  blen 
Stabile  d* apr&s  lea  travaux  do  Kraiclumii  et  Lilley,  et  i’on  obtlent  une 
concordance  raiaonnabie  aveo  lea  experiences  de  Wlllmarth,  Hodgson  et 
d'  autres.  Pour  des  nombres  de  Mach  sieves,  dea  meaurea  existent  grace 
aux  travaux  de  Klstler  et  de  Chen  mala,  Jusqu'a  present,  on  no  dispose 
d’  auoune  thEorle,  4  part  oelle  que  l’on  dolt  k  Phillips,  qul  est  toute- 
foia  apparentEe  au  bruit  Emanant  d’  ecoulements  de  olsalllement  turbulents 
supersonlquea. 

Ce  rapport  passe  en  revue  la  thEorle  des  fluctuations  de  presslon  a 
des  parols  dans  un  Ecnulament  Incompressible  et  montre  comment  la  nature 
dos  fluotuations  da  prssslon  ss  modifle  en  passant  de  1* Ecoulement  &  la 
parol,  On  attire  1* attention  sur  las  interactions  plus  lmportantea 'qul 
donnent  lieu  aux  fluotuations  de  preoalon,  alnal  qu' a  la  rEglon  de  la 
couche  limitrophe,  qul  aausent  prlnoipalement  lea  fluotuations  de  preselon 
sux  psrols, 

On  s  Elargl  les  travaux  pour  englober  les  effets  de  la  compressibility 
On  constate  qu’ une  analyse  slmllslre  1  oelle  de  Phillips  ast  approprlEt, 
blen  que,  oontrsl.'ement  4  oes  dernlers  travaux,  ce  nouveau  traitement  ne 
aolt  pas  llmltE  au  oas  de  nombres  de  Maoh  hypersonlquea  trie  Clevis, 

1/  analyse  se  sort  du  rapport  a_/uT  comma  grand  psramitre  oil  a, 
reprEsente  U  Vitesse  du  son  i  la  parol,  at  ur  represents  la  vi tease  de 
olsalllement,  Cals  est  osrtalnsmant  vrai  pour  une  trie  grande  plage  de 
nombres  de  Maoh,  &  la  condition  qua  la  parol  ne  aolt  pea  soumlse  l  ds 
grandee  vlteases  de  trsnsfert  de  ohaleur.  On  montre  que  les  fluctuations 
de  presslon  sux  psrols  sent  malntonant  la  rEsultat  de  fluotuations  tant 
dans  I'Etat  de  tourbilloiis  que  dans  les  modalltEe  de  son,  Pour  de  grsnds 
nombres  de  Maoh,  oette  dernlEre  contribution  set  sous  forme  d'ondes  Msoh 
de  tourbi lions,  corame  avancE  par  Phillips,  En  faisant  oertalnes  supposi¬ 
tions  conoernsnt  les  Interactions  dominantes,  on  fait  des  Evaluations  de 
la  grandeur  et  du  speotre  des  fluotuations  de  presslon  sux  parols  qul 
rEvilent  des  tendances  slmllslres  aux  mesures  de  Klstler  et  de  Chen, 
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WALL  PRESSURE  FLUCTUATIONS  UNDER  TURBULENT  BOUNDARY  LAYERS 
AT  SUBSONIC  AND  SUPERSONIC  SPEEDS 

O.M.  LI 1 ley* 


1.  INTRODUCTION 

The  worK  of  Heisenberg1,  Obukhov9  end  Batchelor3  hae  shown  that  In  isotropic  turbu- 
lenoe,  the  root  mean  square  fluctuating  pressure  is  given  by 

Jfp1)  -  O.BBfiu1  (1) 

where  u*  is  the  mean  square  value  of  any  fluctuating  component  of  the  turbulent 
velocity,  Uberoi*  has  shown  that  a  relation  of  this  form  exists  in  grid  turbulanos, 
but  the  oonstant  in  Equation  (1)  was  approximately  0,88  over  a  wide  range  of  Reynolds 
numbers, 

The  first  attempt  to  oaloulate  the  fluctuating  pressure  field  in  a  turbulent  shear 
flow  was  made  by  Kraiohnan*  who  found  that  at  the  wall  (subscript  w) 

UlsI  a  2  to  12  (2) 

rw 

where  r„  is  the  local  wall  mean  shear  stress, 

Experiments  in  pipe  flow  by  Willmarth17  showed  that 

^(l>5)  Sr  0.006  iflU 9  (3) 

and  that  the  pressure  field  was  oonveoted  past  the  wall,  at  an  average  speed  of 
0.82  UM  ,  where  U*  is  the  speed  of  the  uniform  flow  external  to  the  boundary  layer. 
However,  in  spite  of  oare  to  reduoe  the  extraneous  noise  in  the  air  supply  leading  to 
the  pipe,  Willmarth  was  unable  to  obtain  acourate  readings  of  the  power  spectral 
density  in  the  lower  frequencies.  Other  investigators  have  obtained  similar  results 
but  only  recently  has  an  attempt  been  made  to  check  the  oonstant  in  Kraiohnan* s 
formula.  (A  more  complete  review  of  the  experimental  work  on  wall  pressure  fluctuations 
will  be  given  in  a  paper  by  Hodgson4  whioh  is  to  be  published  shortly). 


An  additional  problem  noted  by  Willmarth  was  the  correction  necessary  to  allow  for 
the  affenta  n r  r.hfl  finite  size  of  the  oressure  transducer  on  both  the  root  mean  square 
and  power  spectral  density  measurements.  With  this  correction  applied,  most  of  the 
available  measurements  suggest 


1.8  to  3 


(4) 
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over  a  moderate  range  of  Reynolds  numbers  at 

m . 

Tw 

where  a<H)  is  a  slowly  varying  funotion  of 
Reynolds  numbers.  -  " 

Reoent  work  by  Klstler  and  Chen7  has  extended  the  measurements  to  high  Maoh  numbers 
and  their  work  shows  that  a(R)  Increases  progressively  with  Maoh  number,  reaching  a 
value  between  S  and  6  at  a  freestream  Maoh  number  of  5,  at  least  for  the  oase  of  zero 
heat  transfer.  Their  results  suggest  that,  at  a  Maoh  number  of  5,  the  funotion  a(R) 
has  nearly  reached  its  asymptotic  value  for  very  high  Maoh  numbers. 

The  work  of  Kraiohnan  (loc.cit. )  has  been  reviewed  and  extended  by  Lllley'  and  hy 
Llllsy  and  Hodgson9.  The  latter  work  showed  that  the  lower  estimate  of  a(R)  obtained 
by  Kraiohnan  was  more  correot,  and  this  work  also  went  some  way  towards  confirming  that 
the  pressure  fluctuations  in  a  turbulent  shear  flow  are  dominated  by  the  mean  shear. 

The  oaloulated  spectrum  funotion  for  the  wall  pressure  fluctuations  showed  moderate 
agreement  with'  the  measured  spectra  at  high  frequencies,  but  at  lower  frequencies,  the 
oaloulated  fall  was  not  observed  in  the  measurements  made  in  pipes,  wind  tunnels  eto. 

The  corresponding  two-point  pressure  oovarlanoes  showed  marked  differences  between 
longitudinal  and  transverse  separations,  while  the  area  under  the  longitudinal  preesure 
oovarianoe  and  the  related  autocorrelation  was  exaotly  zero.  In  faat,  the  theory 
ehowed,  In  agreement  with  the  work  of  Phillips10,  a  vanishing  surface  Integral  of  the 
two-point  pressure  oovarianoe  taken  over  the  wall.  The  differences  between  theory  and 
the  measurements  of  Wlllmarth  and  others  have  been  investigated  by  Hodgson9,  He  showed 
that  the  ill-defined  strong  negative  loop  in  both  the  measured  longitudinal  pressure 
oovarianoe  and  the  autocorrelation,  and  the  non-vanishing  transverse  pressure  oovarianoe 
at  large  separations,  were  the  result  of  extraneous  disturbances  external  to  the 
boundary  layer.  (The  effeots  of  extraneous  disturbances  were  also  known  to  Wlllmarth 
and  are  also  discussed  at  some  length  in  the  reoent  work  of  Wlllmarth  and  Wooldridge11). 
The  measurements  mads/ by  Hodgson  (loe.clt.),  on  the  wing  of  a  glider  in  flight,  whioh 
were  free  from  extraneous  disturbances,  confirmed  the  relation 


low  Maoh  numbers,  or 

a(R)  (5) 

Reynolds  number  at  surrioientiy  nigh 


i/(p5)  -  2l2  Tw  (6) 

and  showed  the  falling  spectrum  in  the  lower  frequencies,  together  with 

f 

I  P(0;v)dr  =  o  (7) 

\JmOJ 


where  P(0;r)  is  the  autocorrelation  of  the  pressure  at  the  wall.  On  applying  the 
conveoted  hypothesis,  which  is  supported  by  all  the  measurements,  Hodgson  finds  that 
Equation  (7)  Is  equivalent  to 


,.o,£s)df,d£s  = 


0 


(8) 


3 


whioh  Is  an  experimental  confirmation  of  Phillips'  result. 

If  we  return  to  the  problem  of  wall  pressure  fluctuations  at  supersonic  Mach  numbers, 
we  find  that  no  theory  exists,  apart  from  the  work  by  Phillips1*,  on  the  related 
problem  of  sound  generation  by  supersonio  turbulent  shear  layers, 

Phillips  has  shown  that  the  radiated  sound  arises  from  eddy  Maoh  waves  whioh  are 
generated  by  some  wave-numbers  of  the  turbulepoe  in  those  layers  of  the  shear  flow  for 
whioh  the  difference  between  the  mean  velocity  of  the  fluid  outside  and  the  looal  eddy 
convection  velocity  la  greater  than  the  speed  of  sound  outside  the  zone.  Phillips  does 
not  Inolude  the  oase  of  a  wall  shear  flow,  although  clearly  thla  must  present  an 
analogous  problem,  and  Indeed  Phillips  argues  that  his  model  should  be  qualitatively 
oorreot  In  this  case.  However,  measurements  by  Laufer13  of  radiated  sound  from  super¬ 
sonio  turbulent  boundary  layers  are  not  In  good  numerical  agreement  with  Phillips' 
theory,  although  undoubtedly  some  aspects  of  the  phenomenon  desurlUed  by  Phillips, 

Buoh  as  the  production  of  eddy  Maoh  waves,  do  exist  and  have  been  observed  by  many 
workers.  However,  as  Uufer  points  out,  the  experimental  Maoh  numbers  may  not  be  high 
enough  for  Phillips'  asymptotlo  theory  to  be  applicable  in  the  range  of  freestream 
Maoh  numbers  up  to  6.  The  more  general  problem  of  the  sound  radiated  from  shear  flows 
at  supersonio  speeds  has  been  treated  by  Pfowos  Williams14  and  Lighthlll111. 

The  present  paper  seta  out  to  extend  the  theory  of  pressure  fluctuations  in  turbulent 
boundary  layers  In  Incompressible  flow  to  that  at  higher  speeds,  and  to  provide  a  basis 
for  comparison  with  the  measured  results  of  Kistler  and  Ohen  (loc.oit.),  Williams14 
and  Willmarth17. 


2.  INCOMPRESSIBLE  PLOW  T1IE0BV 
2. 1  The  Presiure  Covariance 

It  has  bean  shown  by  Lilley  and  Hodgson  <loc.clt. )  that  the  pressure  at  the  wall  is 
dominated  by  contributions  from  the  turbulence  in  the  'inner'  region  of  the  boundary 
layer,  extending  up  to  about  l.eSj  ,  where  Sj  Is  the  boundary  layer  displacement 
thiokness.  In  this  region,  the  typloal  length  and  velocity  scales  of  the  flow  are 

and  u,  =  /—  respectively  . 

P^r  V  Po 

Measurements  in  this  region  indicate  that  all  except  the  larger  wave  numbers  of  the 
turbulence  are  being  convected  at  a  mean  speed  of  near  0.8  U„,  ,  and  a  theory  of  the 
sub-layer  of  the  'inner'  region,  based  on  this  hypothesis,  is  given  by  Sternberg'". 
Hence  we  might  expeot  that  the  pressure  at  the  wall  ie  also  dominated  by  eddies  having 
this  conveation  speed  of  near  0.8  Um  .  Since  the  correlation  lengths  for  the  wall 
pressure  are  of  the  order  of  a  boundary  layer  thicknese,  It  would  seem  reasonable  to 
neglect  the  rate  of  growth  of  the  boundary  layer  in  calculations  of  the  wall  pressure 
fluctuations.  We  will  assume,  therefore,  that  the  mean  flow  field  is  given  by 
[U  (x  ),0,0]  where  xt  is  measured  in  the  direction  of  the  mainstream  and  xz  Is 
normal  to  the  wall.  If  all  terms  in  the  equations  of  motion  are  made  non-dimensional 
with  respect  to  u T  and  fj.0/p0uT  ,  we  find 


4 


Bu,  dU,_ 

+  Us— *8! 

‘3*!  'dXj 


3(u1u<  -  gpj) 


-iitvs 


where 


<“o 


*1 


PoH^r  . 
Mc 


(8) 


P 


Since  the  Equation  of  continuity  for  the  turbulenoe  la  Buj/Bxj  =  0  ,  we  find,  on 
taking  the  divergence  of  (9),  that  the  Equation  for  the  pressure*  Is 


7»p  s  3dU,  Bu,  -  E^) 

dXj  SiT[  ^XjBxj 


=  A(x,  t)  (10) 

The  two  terms  contributing  to  A(x,t)  ,  which  defines  the  velocity  field,  oan  be 

referred  to  as  the  mean  shear  -  turbulenoe  interaction  (M-T)  and  the  turbulenoe* 
turbulehoe  interaction  (T-T)  respectively. 


The  solution  of  (10)  oan  be  put  in  the  form 

Pftt)  =  '  ~J°  dx'Jdx'dx'  A(x',t)[a0  +  a,] 


Tr, jj* M 


i 

0  Bxi 


(11) 


where  the  surface  integral  is  taken  over  the  wall  at  x'  =  0  .  The  Qreen  functions 
Q0  and  Qj  are  given  respectively  by 

*  It  is  incorreot  to  argue  that  the  second  term  on  the  right  hand  side  of  (10)  is  small  by  oom- 
oarison  with  the  first  term,  However,  integrals  involving  A(s,t)  ore  usually  dominated  by 
the  (M  -  T)  terms  unless  this  contribution  Is  identically  zero. 
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Qo  =  I? 


X 

1-V  ^ 


•  i-» 


Gj  =  |x  -  x'*|*1 


(12) 


where  xf*  ■  (x[,-x',x|)  ie  the  ‘image'  point.  However,  from  the  Equation  of  motion 
(9),  we  see  that  at  the  wall 


3p  _  3aua 

*7  "  3*s 

alnoe  both  U(  and  Uj  vanish  at  the  wall,  and  then  (11)  becomes 


(13) 


a  ‘  dx*jjdx{dx'(G0  +  Qt)  A(x',t) 
®  -a 


(14) 


showing  that  the  pressure  fluctuations  can  be  determined  once  the  velocity  field  Is 
known.  The  covarianoe  between  the  pressure  at  any  point  (x,t)  and  a  random  function 

A* 

q"(x",t),  such  as  the  pressure  or  a  component  of  the  turbulent  velocity  at  (x",t),  Is 

A«>  A* 

therefore  given  by 


1 


p(x,t)q"<x",t)  =  dx' 


47T  | 

1 

27V 


dx{dx'(G0  +  Gj)  A(x')q"(x") 


(16) 


x  a*o 


where  the  bars  indicate  time  means.  Thus  In  order  to  determine  the  pressure  covariance 


p(x)p"(x")  anywhere  In  the  shear  layer  we  need  values  of 


AWSPi  - 

~  'I*  dx'  dXj'dX  j 


(16) 


»  l 

But  p"u,  can  be  determined  from  (15)  by  replacing  q"  by  ui'  so  that 
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—  dx'  dx{dx'(G0  +  Gj) 


dUjSuJuJ  MJuJnJ 

dxj  3x{  SXjfdXj 


*J*0 


(17) 


while  puj'uj 11  la  similarly  obtained  by  replacing  q"  in  (IB)  by  Uj'uJ'  ,  Hence  the 
determination  of  the  pressure  covarianoe  formally  involves  the  evaluation  of  integrals 
of  the  form  of  (15)  and  (17)  over  the  entire  flow  field,  However,  this  oannot  be 
performed  with  any  great  precision  since  the  second-order  veloolty  correlations  are., 
incompletely  known  and  little  is  known  of  the  third-order  and  fourth-order  veloolty 
correlations.  The  mean  square  of  the  pressure  at  the  wall,  p(0)}  ,  and  the  presaure- 

veloolty  covariance,  p(0)u''(x")  are  just  two  particular  results  which  can  be  obtained 

i  ■  A#  . . 

from  the  general  relations  (15)  to  (17). 


On  the  assumption  that  the  <M  -  T)  and  (T  -  T)  terms  are  independent,  Hodgson  (loo, 
clt. )  has  shown,  on  using  the  best  available  data  for  the  mean  and  turbulent  veloolty 
flaw  fields  and  making  extensive  numerical  caloulationB,  that  the  contributions  to 

P(0) a)  ■ 

Tw 

from  the  (M  -  T)  and  (T  -  T)  terms  are  respectively  2,6  and  0.5.  This  shows  that  the 
contribution  of  the  (T  -  T)  terms  to  the  mean  square  of  the  wall  pressure  is  only 
4  per  cent,  and  in  view  of  the  approximate  nature  of  the  oaloulatlons  oan  be  assumed 
negligible.  The  accuracy  of  the  computations  leading  to  the  (T  -  T)  contribution  to 
the  wall  pressure  is  poor  due  to  the  many  assumptions  which  must  be  included  if  a 
numerical  result  is  to  be  obtained.  (The  difference  between  these  calculations  and 
the  measured  values  is  not  considered  to  be  of  major  Importance,  in  view  of  the  fact 
that  on  each  occasion  a  velocity  flow  field  closer  to  the  experimental  one  was  used, 

a  value  of  nearer  to  2.2  was  obtained,) 

If  we  turn  next  to  the  evaluation  of  p(0)u''/^(p(0)2)^/(uj'z)  we  find  that  the 
contribution  from  the  surface  integral  is  negligible  if  x"  »  l  (Note  that  in  our 
notation  x2  a  p0uTXj//z0)t  The  contribution  from  the  (M  -  T)  term  can  be  obtained  on 
choosing  a  suitable  form  for  e  u7u"  (The  contribution  from  the  (T  -  T)  term 

cannot  be  obtained,  even  approximately,  since  values  of  ufu/u"  have  not  been  measured 
except  for  zero  separation).  It  Is  found  that  the  contribution  from  the  (M  -  T)  term 
gives 

P(0)u''  ~  x"F(x")  exp|-/3(x'')(xj,!  +  x"2)*}  (18) 

and  a  comparison  with  the  measurements  of  Willmarth  and  Wooldridge27  is  shown  in 
Figures  1  ana  2.  The  agreement  Is  reasonable,  qualitatively,  except  at  large  separa¬ 
tions,  which  is  nut  surprising  since  the  chosen  form  for  neglected  the 


■  nil-  nil 
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contribution  from  the  big  eddies.  A  Blightly  modified  form  of  (18),  together  with  an 
additional  term  to  represent  the  effeot  of  the  big  eddies,  enables  R„  — «•  to  be 

P  OXj  ' 

represented  in  the  form 


R|)  3x" 


~  f(x 


\ J)  exp 


a(x 


*.*?)}  -  |ax»*)exp{-  J(xf*  +  «;*>■ 


a») 


where  a  is  a  funotlon  of  xj  . 

_ .  <)u!! 

The  evaluation  of  p(0)!  can  now  be  attempted  by  using  R„  — »  defined  by  (IB). 

.  p 
It  oan  be  shown  that  the  contribution  from  the  surface  integral  oan  be  negleoted, 
while  again  it  would  appear  not  unreasonable  to  assume  that  the  (?  -  T)  contribution, 
involving  as  it  does  third-order  oovarlanoes,  oan  be  negleoted.  It  is  found  that 

a/(  PiOp)  has  a  value  of  the  order  of  2,0  if  the  values  of  f(x")  and  a  are  chosen 
appropriately,  The  integral  on  which  this  value  of  <J(p*(0))  is  based  is 


17K 


0(x})dx} 


(20) 


where  0(x9)  Includes  both  the  mean  shear  distribution  and  the  variation  of 
across  the  boundary  layer, 

K  is  the  von  Kftrmdn  constant  used  in  defining  the  mean  velocity  shear.  The  con* 

tributions  to  ^(p][)  arise  from  layers  between  the  wall  and  x}/Sl  =  2.0  ,  and  if  we 
take  a  1610  (corresponding  to  the  value  used  in  one  of  Wooldridge  and  willmarth’ s 
experiments),  this  occurs  at  values  of  x4  up  to  3200,  which  straddles  the  value  of 
xs  =  2550  ,  at  which  Uj/Uo,  =  0.83  .  It  must  not  be  assumed,  however,  that  the  major 

contribution  to  p(0)*  occurs  over  a  small  region  of  the  shear  layer.  In  fact  the 
contributions  are  spread  diffusely  over  a  fairly  large  range  of  x,  as  shown  also  in 
the  experiments  of  Willmarth  and  Wooldridge  and  in  the  analysis  of  Hodgson,  Blightly 
different  conclusions  have  been  suggested  by  Corcos21  recently  from  an  exact  numerical 
integration  of  the  experimental  data  of  Willmarth  and  Wooldridge.  The  result  of  this 

computation  gives  a  value  of  p(0)2)  of  about  1,2  compared  with  the  value  of  2,0 

found  above.  Reasons  for  this  difference  are  not  difficult  to  trace  but  it  is  diffi¬ 
cult  to  assess  which  is  closer  to  the  true  value  of  the  (M  -  T)  contribution  to 

J(p(0)s)  .  Some  objections  to  the  use  of  the  Wooldridge  and  Willmarth  data  at  large 
separations  in  this  computation  can  be  made  on  the  grounds  that: 


(a)  the  experimental  accuracy  is  poor; 

(b)  the  use  of  a  high  pass  filter  In  these  experiments  will  affect  the  values  of 
R  *  at  large  separation?; 


!s .t » tin  *  a-Aid  ikt^  IL  ■  Limit  ‘-1 
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(c)  the  assumed  relation  between  R<  and 

hypothesis' ,  is  not  correct  for  large  separations, 

(o)  is  not  true  for  double  velocity  oovarlanoes  at  large  separations  and  is  likely 
to  be  even  less  true  for  the  velocity-wall  pressure  oovariannn*  measured  by  Wooldridge 
and  Willmarth. 


The  extension  of  this  work  to  include  the  evaluation  of  p2(xg)  involves  many 
difficulties,  although  formally  it  oan  be  obtained  from  (15)  to  (17).  In  view  of  the 
fact  that  3p/3xs  is  nearly  zero  at  the  wall,  together  with  the  result  that  the  major 

contribution  to  p(0)s  arises  from  values  of  x,  up  to  l.sSj,  leads  us  to  suggest 
that  p*  is  nearly  constant,  varying  at  most  by  a  factor  of  2  over  most  of  the  'inner' 
region.  Since .over  most  of  the  constant  stress  layer  u2  1s  of  order  unity  and 

-  puj  ■  -  puj/^u’  is  of  order  4,  we  find  some  confirmation  in  this  suggestion,  *  An 
alternative  suggestion  by  Rentenyik  and  Kovasznay19  is  that  the  fluctuating  pressure 
falls  rapidly  outside  the  'laminar  sub-layer’,  but  this  is  not  in  agreement  with  our 
results.  However,  the  results  in  Section  3  do  ehow  that  the  major  contribution  to  the 
wall  pressure  comes  from  layers  closer  to  the  wall  bb  the  Maoh  number  is  Increased, 
and  ^t  is  probably  this  affect  which  might  have  some  bearing  on  the  results  obtained 
by  RSmenyik  and  Kovasznay. 

1  We  will  now  dleouss  the  pressure-velocity  product  Wt  .  This  term,  whioh  vanishes 
at  the  wall,  is  known  to  play  an  important  role  in  the  energy  transfer  across  the 
boundary  layer,  and  has  roughly  a  constant  value  across  the  entire  ‘inner’  region  of 
the  boundary  layer  outside  the  vlsoous  layer,  Clearly  the  pressure-velocity  covariance 

p(Xj)Uj  must  have  an  essentially  different  form  from  p(0)u"  ,  since,  as  we  have  shown 
above,  p(0)u j(0, Xj, 0)  is  zero,  whereas  pu?(x2)  is  clearly  finite.  In  addition, 

the  surface  integral  of  pu^  taken  over  a  plane  parallel  with  the  wall,  must  vanish 
if  there  is  no  disturbance  outside  the  boundary  layer,  The  results  of  Wooldridge  and 
Willmarth  are  in  agreement  with  this  boundary  condition  for  p(0)u"  . 

The  modifications  to  pUJ  ,  as  the  pressure  measuring  station  is  moved  away  from 
the  wall,  can  be  shown  to  depend  on  contributions  from  the  surface  integral  and  the 
(T  -  T)  term  in  (17),  both  of  which  have  been  shown  to  give  negligible  contributions 

to  p(0)u"  ,  Indeed  if  ^(Xji^)  is  symmetrical  about  =  0  ,  the  contribution 
of  the  (M  -  T)  term  to  pu^  is  zero,  as  noted  by  Corcos20,  and  so  pu^  depends 

*  Ir,  tiia  tut  ui  tha  constant  stress  layer  if  we  nut 


3u. 

Rp  ,  involving  'Taylor’s 


q\  ■  p2v 

and  IpuJI  -  /(p2)v 

then  if  y  -  1 ,  P)  ”  2‘/(p1)  • 
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entirely  on  the  surface  integral  in  (17)  together  with  the  contribution  from  the  third- 
order  velocity  covariance.  Outside  the  viscous  layer,  the  dominant  contribution  to 


PUj  arises  from  the  (T  -  T)  terms.  This  result  is  made  obvious  by  noting  that  if 
structural  similarity  exists  in  the  constant  stress  region  of  the  boundary  layer  and 
outside  the  viscous  layer 


-ss/ 


a,(R)q 
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-Wl 


a,(R)<r 


-;3/s 


and  so 


(21) 


where  q7  is  the  mean  square  or  the  turbulent  velocity,  and  a2(R)  and  aa(R)  are 
slowly  varying  functions  of  Reynolds  number. 


2.2  The  Structure  of  the  Big  Eddies 

The  work  of  Townsend51  and  Grant55  has  suggested  that  the  structure  of  the  big 
eddies  in  a  turbulent  boundary  layer  have  the  form  of  'mixing  Jets'  which  erupt  near 
the  surface  and  spread  into  the  outer  regions  of  the  boundary  layer.-  The  available 
experimental  evidence  in  support  of  this  hypothesis  is  scanty,  although  it  is  consis¬ 
tent  with  Grant’s  own  measurements  of  nine  second-order  velocity  correlations,  and  the 
more  recent  work  of  Wooldridge  and  Wlllmarth  (loc.cit. )  in  which  they  report  extensive 
measurements  of  wall  pressure-velocity  covariances.  Since  the  big  eddies  play  such  a 
vital  role  in  the  determination  of  the  pressure  at  the  wall,  it  is  of  interest  to 
discuss  the  work  of  the  present  author,  in  which  an  attempt  has  been  made  to  put  the 
'mixing  Jet'  hypothesis  on  a  more  quantitative  basis  with  results  in  agreement  with 
the  measurements  both  of  Grant,  and  Wooldridge  and  Wlllmarth.  Only  the  essential 
details  of  this  work  will  be  given  here.  When  the  rateof  turbulent  energy  production 
exceeds  its  equilibrium  level,  it  is  followed  by  Increased  dissipation  and  an  increased 
rate  of  diffusion  of  turbulent  energy  both  towards  and  away  from  the  wall,  The  out¬ 
ward  flux  of  energy  can  be  presented  roughly  by 


U  —  Aq2/2 
dx„ 


where  U  is  the  mixing  jet  velocity  and  Aq5/2  is  the  excess  of  turbulent  energy. 
This  release  of  energy  from  regions  near  the  wall  is  followed  by  an  energy  deficiency, 
and  so  the  outward  'mixing  jot’  must  be  followed  by  a  return  flow  towards  the  nurture. 
Such  a  return  flow  is  also  required  from  considerations  of  continuity.  It  might  be 
expected  that  the  outward  ‘mixing  Jets’  should  have  a  relatively  high  turbulent 
intensity  and  a  relatively  small  scale,  while  the  return  flow  should  have  a  somewhat 
lower  turbulent  intensity  and  a  larger  scale.  An  overall  scale  of  the  ‘mixing  Jets’, 
extending  to  about  a  boundary  layer  thickness,  is  clearly  shown  in  the  measurements 
of  Wooldridge  and  Wlllmarth,  as  well  as  those  of  Grant. 


life),  ;  Mat 
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Since  the  'mixing  Jet'  is  essentially  a  turbulent  flow  away  from  the  surface,  we  ' 
might  expect  the  major  part  of  the  correlation  Rjjfr,,  0, 0)  to  be  the  result  of  the 
relatively  simple  structure  of  the  big  eddy.  The  results  shown  in  Figure  4  confirm 
this.  However  a  fair  fit  with  the  remaining  correlations  measured  by  Grant  could  only 
be  obtained  if,  in  addition  to  the  mixing  Jets,  we  superpose  larger  eddies  rotating  io 
planes  parallel  with  the  wall  and  unoorrelated  with  the  'mixing  Jets*.  The  presence 
of  these  big  eddies  can  be  demonstrated  in  Figure  4,  where  the  separate  contributions 
from  the  ‘mixing  Jets1  and  the  larger  scale  eddies  are  given  for  Rn(r.,  0,  0).  The 
figure  also  shows  the  measure  of  agreement  between  the  model  and  the  measurements  but 
with  a  relatively  free  choice  in  the  values  of  so  many  length  scales  defining  the  three 
part  structure  of  the  big  eddies,  the  agreement  in  many  of  the  examples  shown  is 
probably  fortuitous.  Briefly  the  length  scales  of  the  eddies  have  been  found  to  be  as 
follows: 

(1)  eddies  rotating  in  planes  parallslwith  the  wall  have  a  scale  of  order  8  , 
where  X  Is  the  boundary  layer  thickness,  and  have  a  structure  similar  to 
the  simple  form  suggested  by  Townsend  as  being  representative  of  the  big 
eddies  produced  at  random  in  the  boundary  layer; 

(ii)  the  outward  ‘mixing  Jet'  has  a  scale  of  order  S/io  : 

(ill)  the  scale  of  the  return  flow  is  of  order  8/3  ; 

(lv)  secondary  motions  in  planes  parallel  to  the  wall  accompanying  the  'mixing 
Jet’  have  scales  of  order  S/io  , 

A  diagrammatic  representation  of  the  big  eddy  structure  is  shown  in  Figure  5. 

2.3  The  Pressure  Spectrum 


If  we  define  the  three-dimensional  Fourler-Stleltjes  transforms  of  P(£,t)  and 
A(J(,,t)  respectively  as 


POt-t)  = 

f  i (k,x.+k3x3+6it) 
e  1  1  3  3  deefXgjJi/u) 

(22) 

and 

A«L.t)  = 

f  1  (k.x.+k.x.+iiit) 
e  1  1  3  3  Ay(xi;H,co) 

where  £  «  (kj,k3)  is  the  wave  number  vector  in  the  (Xj,x3)  plane  and  oj  is  the 
frequency,  then  the  equation  for  the  Fourier  coefficient  d w  is 

d£>"  -  kzdw  =  Ay  (23) 

where  primes  denote  differentiation  will)  respect  to  xs  .  If  we  assume  that  the 
disturbance  outside  the  boundary  layer  i3  zero,  d£o(a>)  =  dy(co)  =  0  .  From  the  equation 
of  motion 


dw'(O)  =  dz"(0)  (24) 

where  dz2  is  the  rourier  coefficient  for  u2  . 
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Hence  the  solution  to  (23),  satisfying  the  boundary  condltlon-(24).  is 


4  d£(0)  /kx,  -kxs'| 

dto(x})  =  -  +  6  J 

2 

Ob  ' 


+  j  e‘kydy(y)dy 


-  ^-J  ekydy(y)dy 
and  for  the  Fourier  coefficient  of  the  pressure  at  the  wall 


dw(O)  a  - 


dSJ'(o)  /r 


e"kydy(y)dy 


If  we  assume  that  we  can  neglect  d£!'  (0)  then  the  pressure  spectrum  function  at  the 
wall  is  given  by 

fo  r» 

n<o;k,a>)  =  p  e'!kydy  e"ka  dy(y)dr*(y  +  Z)dz  (27 


y 

where  n(0;^,cu)  =  -ij  P(0;£,t)p(0; £  +  £,t  +  t'je"1^'  £+wt^d^dt' 

877 

J, 

dS(0;k,a))(£>(0;k,cu) 

=  - j.  -  .K -  (28) 

dk^dkjdtu 

and  £s  (Tj, rg)  and  the  asterisk  denotes  the  complex  conjugate,  The  mean  square  of 
the  pressure  at  the  wall  is 


■  Jlf' 


P*(0)  =  J  Ij da)  ri(0;£,a)) 
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Since  all  wave-numbers,  apart  from  the  highest,  are  converted  at  a  constant  Bpeed, 
U0  ,  the  frequency  speotrum  can  be  obtained,  relative  to  co-ordinates  fixed  in  the 
wall,1  from  the  integrated  wave  number  speotrum  funotion.  Thus 

ll(0;kt)  = 

where  11(0;}$)  = 

A  result  of  the  vanishing  surface  integral  of  the  pressure  covarianoe  leads  to 

n<o;  0)  =  o  (31) 

In  these  results  k  and  co  are  dimensionless,  and  are  given  by  t)ie  relations 

kS, 

k  =  - — — 

Pourh^O 

0)8.  /u 

and  tn  =  - L— ! —  ,  (32) 

PourSi^o 

If  A(j, t)  in  (10)  is  Independent  of  Reynolds  number,  then  n(0;k,o)  should  be  a 
universal  function  of  k  and  u>.  (Hence  spectra  plotted  as  functions  of  aJ&l/Um  , 
say,  will  not  be  universal  functions). 

If  we  now  replace  d7(x2;}{,,co)  by  its  (M  -  T)  term  only 

and  substitute  into  Equation  (27) 


n(0;)$/»)  = 


where  Is  the  four-dimensional  spectrum  function,  £  =  (x, ,«,,«,)  ,  *,  =  k,  , 

k3  =  k3  and  r  =  dl^/dXj  . 


«U  a 

-21k.  —i  dz .(XjiJt.o) 

1  dXj  22 


(33) 


4k  * 

l  e‘2kyr(y)dy 
kz 

Jo 


P 


*22(y''£‘0J)d/<2 


•  e~kzei/<!Zr(y  +  z)dz 

J-y 


(34) 


r* 


ll(0;k)dka  =  nfo:  -  — 


(30) 


i: 


n(0;fc,co)d£u 


If  we  assume  that 


*  «  3  \  112/ 


UJ 


^u*mju|(y  +  z)  ^8/*  u0a  +  £***) 

where  Z2k2  ■  Z2k3  +  Z|k| 
as  shown  in  the  Appendix,  and 

.  .  zj u]  =  h(0)e*^  ;  lt  =  y  ;  Zj  =  Z3  =  i 


(35) 


we  find 


n<o;j<,£o)  .=• 


[e'1  +  e2El(-2)  +  e~26l(2)  -  e~2il(l)] 
tr3/2(k  +  /S)*U 


■  3kJh!(0)Z8 


(!2k2+i2a.2/uJ) 


(36) 


where  the  tern  in  square  brackets  Is  equal  to  0.42  approximately.  (The  contribution 
to  n<0;fc,co)  from  the  sub-layer  can  be  shown  to  be  negligible  and  since  the  major 
contribution  is  found  to  come  from  layers  between  (4  to  2S\  ,  we  are  Justified  in 
putting  h(0)  as  finite).  The  Integration  over  all  frequencies  can  be  performed  and 
the  wave  number  spectrum  is  then  given  by 


n(0;k> 


1.68k2ZV(0)e~i8|t2 
7r(k  +  /3) 8 


(37) 


and 


Ps(0) 


1.08hz(O)  Z* 


k3  e_)<  dk 
(k  +  PI  )2 


vo 


(38) 


An  analysis  of  the  contributions  to  17(0;^)  from  different  values  of  y  shows 
that  the  larger  contributions  come  from  layers  of  order  S,  from  the  wall,  at  least 
for  the  energy  containing  eddies.  If  we  then  choose  values  for  h(0)  and  /?  to  give 
a  good  fit  to  h(y)  near  y  =  ,  we  find 

him  =  R.4/S.  :  l  =  8  :  BS.  -  1 
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The  Integral*  can  be  evaluated  In  terms  of 


J(x) 


k  +  x 


Jo 


which  Is  tabulated  by  Qoodwln  and  Staton23. 
We  find,  finally,  that 


K  =  2-2 

which  shows  good  agreement  with  the  experimental  results. 


(39) 


3.  COMPRESSIBLE  FLOW  THEORY 


3. 1  The  Preasuro  Disturbance  Equation  (zero  heat  transfer) 


since  the  experimental  results  of  the  wall  pressure  fluctuations  by  Klstler  and 
Chen  (loc.cit.)  at  high  Mach  numbers  show  only  a  relatively  small  divergence  from  the 

linear  relation  between  JpJ  and  tw  we  might  expect  that  the  dominant  terms  con¬ 
tributing  to  the  wall  pressure  in  Incompressible  flow  also  play  a  dominant  role  In 
compressible  flow.  Thus,  if  we  neglect  the  diffusive  terms  and  write  our  equations  in 
dimensionless  form,  we  find,  following  Phillips  (loc.cit,),  If  the  mean  values  of  the 
density  and  viscosity  at  the  wall  are  constant,  that 


a! 

U1  +  +  .  !l 

ao  DtJ  a2  dx2  Bx2 


(40) 


where 


& 


W l. 

Mw 


tu^°.  . 

^w 


* 


-k 


(It  *  Ely 


dk 


<»hi  -  2/w  d/-)  *  3 (pi)2  j(pl)  *  (4<r  J’ipl) 
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A (l,  *■) 


\  Bxj  av 


(41) 


D  3  3 

—  =  ^-  +  U.  ^ —  , 
Dt  3t  1  Bxj 


where  Ti t  =  U,(x8)  . 

The  speed  of  sound,  a  ,  la  also  a  function  of  xs  only, 


As  Phillips  remarks,  the  left  hand  side  of  this  equation  neglects  convection  and 
scattering  of  the  sound  by  the  turbulence  and  by  fluctuations  In  the  speed  of  sound. 

The  diffusive  terms  can  be  added  to  A(£,t)  and  hence,  by  using  (40)  as  the  baslo 
equation  defining  pressure  fluctuations  in  a  turbulent  aompresslble  flow,  there  Is 
little  loss  of  generality.  We  see  that  the  effect  of  lnoludlns  fluctuations  In  density 
(sound  waves)  has  given  rise  to  additional  terms  as  compared  with  the  equation  In 
Incompressible  flow,  Since  diffusive  effects,  however,  have  been,  negleoted,  we  find 
that  the  pressure  fluctuations  are  the  result  of  the  fluctuating  vortioity  and  sound 
modes,  where  In  general  the  vortioity  mode  is  the  larger.  If  this  were  not  so,  it 
would  imply  that  the  mean  properties  of  the  turbulence  in  a  compressible  shear  flow 
could  not  be  derived  from  a  transformation  of  the  results  in  incompressible  flow. 

But  both  Morkovln84  and  Coles81  have  shown  that  this  scaling  up  of  the  Incompressible 
data  gives  fair  agreement  with  the  limited  measurements  made  in  supersonic  flows.  If 
the  sound  mode  could  be  Ignored,  we  could  write  (40)  In  the  form 


and  its  solution  would  then  follow  on  the  lines  given  above  for  incompressible  flow. 
However,  at  this  stage,  we  will  continue  to  retain  (40)  In  full  and  investigate  more 
fully  the  terms  giving  rise  to  the  sound  mode. 


The  equation  for  the  Fourier  coefficient  dec  ,  as  defined  in  (22),  is,  If 


=  kj  +  kj 


a8 

d  ln-j 

dcS"  +  - -  iuf  - 

dx. 


k  ~  a*-  <£lJ+  uiki)! 


dSS  =  —  Ay 

P w 


(42) 


since  a!/a£  =  pm/p  for  constant  mean  pressure  over  the  entire  shear  layer,  and  primes 
ueiiula  uiftextmLi.ttLi.uii  nriLit  iowyoCt  to  Xg  . 

If,  following  Phillips  (loo.cit, ),  the  first  derivative  Is  eliminated  by  the  use  of 
the  new  dependent  variable 


S'  Arf 

. —  da) 
a* 


Hu»n 


16 


v  -  [ks  -  -> + uk,)2  +  -lc  =  fA_y 
L  a  aJ  yv./ 


with  £(0;}^,to)  =  dto(0 ;&/*))  .  This  equation  should  be  oompared  with  its  corresponding 
equation  In  Incompressible  flow  (23),  which  Is  obtained  from  (43),  when  a  =  aw  -  oo  , 


For  the  case  of  zero  heat  transfer,  we  have 


&S  = 

a;  =  0  ; 

and 


It  can  therefore  be  seen  that  for  snijill  wave  numbers,  a"/a  is  a  not  unimportant  term 
in  the  left  hand  side  of  (43),  However,  except  at  the  wall,  it  is  small  compared  with 
the  source  terms  on  the  right  hand  side  and  so  in  general  can  be  neglected. 

Now  in  his  analysis  Phillips  chose  non-dimensional  co-ordinates  such  that  the  width 
of  the  shear  layer  was  unity  and  the  Mach  number  of  the  external  flow  was  very  high. 

He  found,  finally,  the  radiation  of  sound  from  the  shear  layer  by  a  solution  which 

neglected  terms  of  order  1/Ma,  . 

In  our  problem  we  have  chosen  boundary  layer  co-ordinates  such  that  Uj  -•  Um  as 

x g  “*  co  ,  We  cannot  strictly  estimate  the  radiation  of  sound  from  the  boundary  layer 

since  we  cannot  enter  the  far  field  outside  the  boundary  layer.  However  we  can 
estimate  the  disturbance  in  the  outer  region  of  the  layer  if  we  find  the  solution  to 
(43)  satisfying  appropriate  boundary  conditions.  By  neglecting  diffusive  effects,  our 
boundary  conditions  are 


a"  y-i  u  a*  . 

t  ■  -—iA0'0;*?”!, 


d2)'(0)  =  £'(0)  -  0 
and  either  dcu(co)  =  £(co)  =  0 

or  £  ~  e-lXx2  Xj  on 


(45) 

(46) 


where  we  differentiate  between  the  cases  of  zero  disturbance  outside  the  boundary  layer 
and  that  of  outward  propagating  waves, 


In  our  problem  it  is  too  restrictive  to  find  only  a  solution  for  large  values  of 
and  it  is  desirable  to  choose  some  other  parameter  which  defines  the  flow.  In 
(43)  wo  notod  the  oxistcnco  of  the  term  u2/a2  and  wc  find  that  for  zero  heat  transfer 
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for  all  Mffl*  .  It  follows  that  a  solution  to  (43)  Is  required  for  large  values  of  the 
parameter  A.  =  a,/!!,.  .  (This  is  true  also  for  most  oases  with  heat  transfer). 

Let  us  now  define  the  new  independent  variable 


Ur 


Then  with 


k8aj|  /du  +  U.k.V  5.v  u? 

q(5°  =  +1^ 


(47) 


where  a  =  a/a„  ,  we  find  (43)  becomes 


d*C 

77  -  AJQ(y)5  =  H(y) 
dyz 


(48) 


where 


H(y) 


dy 


and  the  term  in  In  H(y)  is  absorbed  finally  In  the  velocity  derivatives  and  wave 
number  terms. 


If  we  now  follow  the  arguments  used  in  the  incompressible  flow  theory,  we  see  that 
H(y)  will  be  negligible  over  the  outer  three  quarters  of  the  boundary  layer.  Also, 
in  this  same  region,  5yy/5  will  be  small,  owing  to  the  small  gradients  in  the  mean 
velocity.  Thus  a  good  approximation  to  (43)  in  the  outer  region  of  the  boundary  layer 
will  be  to  replace  Uj  by  U„  ,  a  by  aB/aw  and  then 


d_!i  *8 

dy8 


s=s 


K  &U  /  ft-  «A— 

— r  - [u  -1  + k ,h„  -1 

*  \  a  1  11 


(49) 


But  if  the  convection  speed  of  the  turbulence  is  Ue  ,  the  frequency  u  is  given  by 


® 


(80) 


•If,  following  Coles, 


for  the  case  of  zero 


we  put 


heat  transfer,  where 


a8  2(1  -  Kc?i) 


m£(1  -  Kcfl)J 
K  y  153  ,  we  see  that 
as  Mj,  -•  «  . 


.tail* III 4. 
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and  (49)  becomes 


having  solutions  of  the  form 


4  exp(±  qjy)  for  qt  >  0 

£  ~  exp<±i\|q,|ty)  for  qt  <  0 


Q(  =  q(“)  = 


The  interpretation  of  the  Rsnond1  t.yps  of  solution  was  given  by  Phillips,  who  showed 
that  it  is  equivalent  to  toward  and  outward  propagating  eddy  Mach  waves  respectively 
having  wave  numbers  such  that,  with  k}  =  k  cost?  , 


Hence  eddy  Mach  waves  are  generated  over  wave  numbers  in  the  plane  for  which 
77  -  0n  <  9  <  v  +  0n  and  -  dg  <  S  <  0n 


9m  =  cos"  ^  (• - - -  (55) 

™  V1  "  V*V*W. 

and  when  U0/Ua  =  0. 8  ,  outward  radiating  Mach,  waves  occur  when  M  exceeds  1.25  , 

The  region  near  |'Jt|  =  0  is  excluded  since  eddies  beyond  a  certain  size  do  not  exist. 

Although  these  considerations  give  us  some  idea  of  the  solution  to  (48),  they 
really  only  help  us  define  conditions  which  have  to  be  satisfied  near  the  outer  edge 
of  the  boundary  layer.  It  is  shown,  therefore,  that  the  first  boundary  condition  of 
(48)  is  applicable  for  wave  numbers  defined  by 

77  -  d„  <  6  <  and  77  +  £L  <  0  <  2v  -  8m 
m  m  m  m 

while  the  second  boundary  condition  applies  for  all  other  values  of  6  . 

In  regions  closer  to  the  wall,  q(y)  changes  sign  when* 


-  U,)‘ 


(rr‘X".J-' ♦ <  t) 

k2(U0  -  Uj)z 


*  The  terms  involving  a"/a  are  Included  for  convenience  only. 

As  already  stated.  It  is  preferable  put  them  nn  the  right  hand  aide  of  (4fl)  and 
to  regard  them  as  source  terms. 
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and  this  occurs  at  the  wall  when 


Hence  for  y-  1.4  and  Ug/Uj,,  =  0.8  to  0.7  ,  this  limiting  value  of  ranges  from 
1.5  to  1.8  approximately,  provided  that  t)ie  bracket  term  Is  0(1).  At  low  supersonic 
Mach  numbers  there  still  exists  an  extensive  region  of  the  boundary  layer  where  the 
flow  Is  subsonic  relative  to  the  speed  of  sound  at  the  wall,  even  though  the  local 
flow  Is  partly  supersonic,  and  eddy  Mach  waves  do  not  exist.  At  higher  Mach  numbers, 
as  shown  by  Kistler  and  Chen,  the  conveotion  velocity  is  supersonic  relative  to  the 
speed  of  sound  at  the  wall  and  at  each  region  distance  y  from  the  wall,  there  Is  a 
range  of  wave  numbers  from  which  eddy  Mach  waves  are  generated.  The  remaining  wave 
numbers  in  the  turbulence  produce  disturbances  near  the  wall  of  the  exponential  type, 
Just  as  in  incompressible  flow.  The  condition  for  the  generation  of  eddy  Mach  waves  Is 


Even  so  except  at  very  high  Mach  numbers,  there  is  a  large  range  of  wave  numbers  for 
which  the  effective  speed  of  the  disturbances  is  subsonic  with  respect  to  the  wall. 

3.2  The  Solution  of  the  Pressure  Equation 

We  have  shown  above  that  the  pressure  disturbance  equation  can  be  written  in  the 
form 

d!£ 

— -j-  -  A-*<i(y)£  =  H(y)  (48) 

dy! 

where  q(y)  and  its  derivatives  are  continuous  functions  of  y.  This  equation  has  a 
transition  point  at  y  =  y0  ,  where  q(y0)  =  0  .  If  y  =  Y  when  o>  +  U,kt  =  0  where 
cu  =  -Ujkj  ,  then  q(y)  will  be  positive  in  the  range  y0  <  y  <  Y  ,  provided  that 
Eyy/I\8  can  be  neglected.  If  a  second  transition  point  occurs  at  y  =  where 
q(yj)  =  0  and  yj  >  y0  ,  then  q(y)  will  be  positive  in  the  range  y0  <  y  <  yt  ,  but 
q(y)  will  be  negative  for  0  <  y  <  y0  and  y  >  yl  . 


Case  I  q(y)  >0  0  $  y  $  ys 

Since  \  is  a  large  parameter,  we  can  find  an  asymptotic  solution  to  (48).  On 
neglecting  terms  of  0(lA)  and  Inserting  the  boundary  condition  d£/dy  =  0  at 
y  =  0  we  find 
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where  A'  is  a  constant  given  by 


{(0)  =  d£o(0)  = 


but  can  only  be  determined  when  some  other  boundary  condition  is  inserted,  for 
instance  at  y  =  y2  .  If  we  argue,  however,  that  the  region  surrounding  y  =  Y 
provides  the  dominant  contribution  to  the  disturbance  at  the  wall  and  for  y  >  y2  , 
£(y)  «  1(0)  we  find  that 


1  H(y'l 

Xq(0^  d(y') 


exp^-\  Jy  q^dy"jdy' 


v  V 

which  reduces  to  (26)  in  incompressible  flow  when  =  — —  and  y0  =  oo  . 

U* 

Cose  II  q  <  0  yt  <  y  <  a> 

If  a  second  transition  point  occurs  at  y  =  yi  ,  then  q(yj)  =0  and  q  is 
negative  in  the  region  <  y  ^  oo  , 


Let  us  put 


q*(y)  =  -q(y) 


/co  +  U jk j\  2  ayy  u£ 

V  is-/  _  “ul  if  a* 


then  a  solution  is  required  of 


— j  +  ^V(y)£  =  H(y) 


where  q*(y)  >0  in’  y1  <  y  <  oo  and  q*(yj)  =  0 


Near  y  -  we  have 


where  u* '(!/,)  >0 


q*  ~  q*f(y,Ky  -  yt) 


If  we  introduce  the  new  dependent  variable  fi(f)  where 
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and 


*  =  q‘* 


(66) 


then  from  (63) 


d*n 

— r  t  X~gVl 
d£* 


with 


a 

(f'>* 


+  small  terms 


f  = 


q^dy' 


i/3 


The  range  of  £  is  0  S  £  <  00  and  we  take 


d£  /q*yS 

<«y  =  \  f) 


If  we  assume  that  only  outgoing  Mach  waves  exist  for  y  »  y  ,  then 


(67) 


(68) 


(69) 


£(y)  ~  exp(-i\  /q*dy' 


(70) 


The  solution  of  the  homogeneous  part  of  (67)  can  be  written  down  In  terms  of  the 
funotlons 

M*,/3  J*i(M*> 

where  /x*  =  §A.<f3/2  and  0  £  fjf  <  co  and  the  solution  for  £(y)  ,  which  can  be  con¬ 
tinued  analytically  into  the  region  y  (  yj  ,  follows  apart  from  one  undetermined 
constant. 

Case  III  q  >  o  y0  <  y  <  y,  * 

As  stated  above,  the  solution  obtained  In  the  region  yj  <  y  <  oo  can  be  continued 
analytically  Into  the  region  y0  <  y  <  y;  and  therefore  involves  terms  of  the  form 

M1/31 i(AO  and  M1/3Kl(M) 

3  3 

where  n  =  §\t3/2 

~3  fy  1  12/3 

and  t  =  _  I  i/q"dy' 

L2  Jy  J 


The  solution  for  £(y)  Is  therefore  given  in  this  region  also,  apart  from  one 
undetermined  constant,  the  same  constant  as  in  Case  II  above.  This  solution  cannot, 
however,  be  extended  to  values  of  y  near  y0  where  q(y0)  =  0  .  A  solution  can, 
however,  be  found  around  y  =  y0  involving  the  functions 
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where 

and 


Vil3IUV)  and  Vil3^l(V) 

3  3 

V  =  §Xs3/2 


The  solution  contains  two  undetermined  constants. 


Case  IV  q<0  0<y<y0 

The  solution  obtained  in  111  can  be  continued  analytically  into  the  region 
0  <  y  <  y0  ,  and  hence  it  is  given  apart  from  the  two  constants.  Boundary  conditions 
at  the  wall  provide  one  relation  and  two  more  relations  are  obtained  by  matching  the 
solutions  obtained  in  III  at  some  convenient  value  of  y  In  y0  <  y  <  Sl  .  The  three 
unknown  constants  can  now  be  evaluated  and  the  solution  for  £(y)  throughout  the 
entire  boundary  layer  has  been  obtained,  and  in  particular  the  value  at  the  wall. 

(Pull  details  of  this  solution  are  given  in  a  separate  paper  which  will  be  published 
shortly). 

First  of  all  let  us  find  the  changes  that  result  in  d25(0)  as  a  result,  of  com¬ 
pressibility  effects  at  low  supersonic  and  high  subsonic  Mach  numbers  and  these  can  be 
demonstrated  on  evaluation  of  (81), 

If  we  assume,  as  in  incompressible  flow,  that  the  dominant  contribution  to  H(y) 
arises  from  the  (M  -  T)  term  then 


H(y) 


(71) 


where 


a2 

d(y)  ^  k2  — 


(co  +  Ujkj)2 

a2 


We  now  require  values  of  dVl/dx2  and  u2  as  functions  of  the  freestream  Mach 
number,  and  these  can  be  obtained  from  their  so-called  equivalent  incompressible 
counterparts  using  Coles’  transformation  formulae. 


If  quantities  with  a  bar  represent  incompressible  values 


(73) 


where  pe  is  the  sub-layer  viscosity  which,  for  zero  heat  transfer  and  y  =  1.4  ,  is 
evaluated  at  the  temperature  Tg  given  by 


(1  +  o.um£) 

(1  +  0.2M|) 


(74) 


Also 


(75) 


where  0t  is  the  non-dimensional  mean  velocity  in  incompressible  flow  and  Uj 
is  its  corresponding  value  in  compressible  flow.  If  further  we  assume  that  the 
relation  between  T  and  Uj  is  given  by  the  Crocco  energy  integral,  we  can  perform 
the  integration  in  (73)  and  so  find  the  relation  between  x2  and  S2  in  the  form 


y  - 

1  M2 

x  -  5  fie 

2 

.*}] 

2' 

'  *K 

1+I 

71* 

'J«t 

VJ 

where  K  is  the  von  K&m&n  constant.  In  addition  it  is  found  that 


Pa  dUj 

'  a r2 


dUt 

ixT 


r  y  - 1 

M2 

3  /  2 

-2  2 

Moo 

n# 

1  u, 

i  +  y- 

1  2 

L  2 

Mco 

_ 

(77) 


for  the  case  of  zero  heat  transfer. 

If  we  assume  that  uj  changes  with  compressibility  in  the  same  way  as  the  mean 
flow* 

uf(xs)^  =  if(V  (78) 

Pa 


*  This  relation  between  u|  and  tT2  differs  from  that  used  by  Morkovin  (loc.oit.) 
but  neither  relation  is  in  good  agreement  with  the  available  experimental  data,  All 
that  can  be  said  of  (79)  is  that  it  qualitatively  has  the  right  trend  with  Increase 
in  Mach  number. 


and  a  similar  relation  is  assumed  to  exist  between  the  Fourier  coefficients  such  that 


dz^Xj)  =  d2j<*j)  ^  (19) 

After  some  reduotlon  we  find,  on  substituting  (76),  (77)  and  (7B)  into  (72),  that 


doi(0;jt,«i)  = 


2lkl  —  JR  , 

1  dll ^  dz2  a 

kX(0)^  dX  ~ 


J. 


Sd*2 


(80) 


where 


vt  -  uS  _  ts 


K*x<x.y-*  =  k*  -  u 


®  +  ki  um/J 


u®  ,  ,  2L^i  u. 


1+  "2  M“/1  ufs? 


and 


fx* 


kf(i-)  =  k 


X(X')  15 


Jo 


M1  4^4 


The  spectrum  function  is  therefore  given  by 


ri(o;k,cu)  = 


r(y)  e*2k£W> 


X(y)^J 


p- dy  .  $j2(y,z;)i,K2,cu)d/<2 


Itfz  -k  Cf  ly*zl-f  fy) ) 

r(y  +  z)e  e 

t  /T(y  +  z)\  ^ 

X(y  +  z)^  ' 


i-y 


dz 


(81) 


where  the  upper  limit  in  (80)  has  been  replaced  by  infinity.  We  see  that  (81)  is 
identical  with  (34),  apart  from  the  term  in  •  T/Tw  •  kX(y)  in  place  of  k 
and  kf(y)  in  place  of  ky  . 
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A  rough  approximation  to  (81)  is  found  by  replacing  f(y)  with  fy  and  X(y) 
with  7  ,  where  1  and  X  are  independent  of  y  .  We  then  find  that  the  integrals 
have  the  same  form  aa  in  Incompressible  flow  and  have  the  values  obtained  previously. 
Hence,  on  reference  to  (35)  and  (36),  we  find 


1.68kjil,hs(o/—  ) 

^  e‘!V  (82) 

where  the  integration  over  all  frequencies  has  been  performed  on  the  assumption  that, 
except  at  very  small  wave  numbers, 

X(y)  *  x  *  X(xs  at  ut  =  U0)  =  1  ;  h(y)  =  =  • 

The  integration  over  all  wave  numbers  and  k3  can  next  be  obtained  and  leads  to 

Pa(0)  =  i.es^  h2(0>  (83) 


where 


1  k3e*k2dk 


(84) 


which  can  be  evaluated  in  terms  of  the  function  tabulated  by  floodwin  and  Staton 
(loc.cit, ). 


All  that  remains  is  to  choose  suitable  values  for  h(0)  ,  fil  and  f  .  If  we 
follow  the  arguments  used  in  the  incompressible  flow  analysis,  we  must  choose  h(0) 
such  that  h(y)  is  a  good  approximation  to  its  corresponding  value  in  compressible 
flow  near  x2  =  . 


It  follows  that  some  adjustment  to 
1  . 


h(0) 


/SSj  with  Mach  number  is  necessary  and  we  put 


(85) 


4 


with 


"MiHfrSfli  T.  iwi-wJH  ...Hi  j  ^  !ipla,i|^ 
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approximately. 

If  we  use  the  asymptotic  expansion  for  the  function  J(x)  given  by  Goodwin  and 
Staton  we  find 

/p{(0)  "*  6.3  as  Mm  -  co  , 

Figure  6  shows  results  for  /p*(0)  evaluated  from  (83),  together  with  (84)  to  (87), 
and  the  experimental  results  of  Klstler  and  Chen  (loc.clt. )  and  Willmarth,  Hodgson 
and  Mull  and  Algranti26  at  low  speeds.  As  already  stated,  these  theoretical  results 
should  only  be  applicable  for  low  supersonic  Mach  numbers  but  it  Is  here  that  we  have 
the  greatest  divergence  with  the  experimental  results.  There  does  not  appear  to  be 
any  Justification  for  reducing  the  experimental  results  of  Klstler  and  Chen,  by  a 
constant  factor,  but  if  this  is  done,  they  fit  both ’the  low  speed  data  of  other 
workers  and  the  theoretical  curve  throughout  the  entire  range  of  Mach  number. 

However  on  the  assumption  that  the  difference  between  the  present  theory  and 
experiment  is  real,  clearly  we  must  find  what  is  wrong  with  regard  to  the  theory.  It 
1b  difficult  to  see  how  it  is  wrong  at  low  supersonic  Mach  numbers  and  high  subsonic 
Mach  numbers,  where  no  eddy  Mach  waves  can  exist  and  where  the  compressibility  effects 
on  the  mean  flow  field  and  the  turbulence  are  known  to  be  small. 

If  we  ignore  the  experimental  point  of  Klstler  and  Chen  at  Mm  =  0. 6  ,  and  assume 
that  at  some  value  of  above  1.25  a  large  increase  in  pressure  level  occurs  as  a 
result  of  the  generation  of  eddy  Mach  waves,  then  it  is  surprising  that  an  ever 
increasing  divergence  between  our  theory  and  experiment  does  not  exist  as  Mm  is 
increased.  The  fact  that  both  sets  of  results  appear  to  have  the  same  asymptotic 
behaviour  at  high  Mach  numbers  seems  to  suggest  that  eddy  Mach  waves  do  not  contribute 
greatly  to  the  wall  pressure  fluctuations.  The  significance  of  this  will  be  explored 
In  the  next  section. 
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Before  closing  this  section  we  note  that  the  changes  In  the  frequency  spectrum  with 
increase  in  will  not  be  large  and  the  peak  should  occur  at  near  =  0.3  , 

Its  value  In  incompressible  flow.  In  fact  the  results  of  Kistler  and  Chen  are  In  good 
agreement  with  the  low  speed  resultB  of  Willmarth,  Also,  the  main  effects  of  Mach 

number  on  appear  to  be  a  reduction  due  to  the  decrease  In  c{  with  Increase  In 

Mo,  ,  and  an  Increase  due  to  the  shift  of  the  dominant  source  region  nearer  to  the  wall 
with  increase  In  , 


3.3  The  Pressure  Equation  at  High  Mach  Numbers 


We  have  shown  above  that  the  solution  given  in  (61)  is  restricted  to  the  oase 
q (y>  >  0  .  If,  therefore,  a  transition  point  exists  between  the  wall  and  the  station 
where  Uj  =  u0  ,  ws  must  turn  to  the  solution  outlined  in  Case  III  above.  The  solution 
follows  the  approach  used  by  Phillips,  although  we  find  it  necessary  to  modify  that 
treatment  when  applied  to  our  problem.  We  will,  however,  still  not  consider  the  full 
solution  which  must  include  radiation  outwaVds, frpm  the  boundary  layer. 

The  solution  to  (48)  in  the  region  around  y  =  y„  where  q(y0)  =  0  is  found  to  be 


where 


i(y) 


1 

rflMv)  o-(y)  + 

3 


T)3Kier))/3(y) 

3 


(88) 


and 


Q(y0)  =  0,  7)=(§)^s3/2 


«(y) 


.i 

<^>(s)773Ki  (rj)ds 

3 


(89) 


and 


£<y)  = 


1 

0(S)T)3Ii(7))dS 

3 


0 


4>(  s) 


(90) 


and  sines  <j>(a)  contains  £(y)  (88)  is  an  integral  equation  for  £(y)  . 
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where 


/3*(0)  =  B 


(1 


<5^(  s* )  rf  3J|(T7*)da* 
3 


and 


(96) 


(0)  =  A - j 


fix.3 


4>( 8*)  v*3(^x(v*)  +  J-iO?*)^ds* 


It  Is  convenient  to  approximate  to  these  integrals  by  assuming  :‘4>(BVY ='"rfbiWt'aftt  =" 
<fi(Q)  ;  then 


a*(0)  ~  A 


A 


H(y„) 

1  A'(y„) 


(97) 


and 


/3*( 0)  ~  B 


19;  .  »(y0) 
J 


(98) 


noting  that  q(y0)  -  o  ,  and  that  for  the  value  of  7)J  given  by  1  =  J|(7)J)/J.| (?)$) 
we  must  put  a*(0)  =  0  in  order  to  satisfy  (94),  and  then  /3*(0)  is  not  determined. 
But  clearly  this  only  arises  because  the  approximation  given  by  (94)  is  inadequate, 
l.e.  higher  order  terns  must  be  retained,  and  therefore  in  what  follows  7$  must  be 
less  than  unity, 

Now  the  value  of  q*  .  Is  found  from  (47)  to  be 

q*  =  cJ  -  kV  +?— i 
2 


o  <  e  <  en 


• - vnn. 


(Tl 


q;  >  0  for 


where 
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which  shows  that  at  9  =  the  speed  along  the  normal  to  the  wave  fronts  Is  roughly 
equal  to  the  speed  of  sound  of  the  gas  at  the  wall.  Thus  the  position  of  the  transi¬ 
tion  layer  at  which  y  =  y0  ,  q(y0)  =  0  ,  changes  with  frequency  and  coincides  with’1 
the  wall  when  q*  =  0  .  In  the  range  0  <  6  <  6n  ,  qB  <  0  ,  while  in  the  range 
9n  <  0  <  77/2  .  qw  >  O'.  For  the  latter  region  we  can  make  use  of  the  solution 
obtained  in  (60) 

2A  ^  77 

or  8  dc3(0)  =  -y  6a  <  9  <  -  (60) 

2 

The  value  of  £w  when  qj  =  0 


d2)(0)  = 


Similarly,  from  (95) 


do>(0)  = 


So  far  we  have  applied  only  the  boundary  condition  that 'the  normal  pressure 
gradient  vanishes  at  the  wall  when  diffusive  effects  can  be  neglected.  However  one 
further  condition  is  required  in  order  to  determine  the  constants  in  the  formulae  for 
4(0)  .  This  further  boundary  condition  results  from  the  disturbance  level  near  or 
beyond  the  outer  edge  of  the  boundary  layer.  If  we  assume  that  the  bulk  of  the  pressure 
disturbance  near  the  wall  arises  from  the  region  around  y  =  Y  ,  it  would  seem  not 
unreasonable  to  assume  that  the  level  of  the  pressure  fluctuations  near  y  =  yt  , 

Miere  yt  >  Y  and  q(yt)  =  0  ,  must  be  very  much  smaller  than  at  the  wall,  in  spite 
of  the  fact  that  radiation  outwards  is  taking  place,  As  already  stated,  the  proper, 
boundary  conditions  at  the  edge  of  the  boundary  layer  can  only  be  applied  to  the  region 
y;  <  y  <  co  ,  where,  for  sufficiently  high  Mach  numbers,  there  will  exist  a  range  of 
wave  numbers  for  which  q  will  be  negative.  However,  for  our  purpose,  it  will  be 
sufficient  to  assume  that  y0  and  yt  are  well  separated  so  that  we  can  put 

4  =  0  for  y  »  Y  . 

(Of  course  the  ,radlation  outwards  can  only  be  determined  by  proper  matching  of  the 
solutions  around  y  =  y0  and  y  ~  yl  as  previously  discussed),  we  therefore  find 
that  a(y)  =0  for  y  »  Y  or  from  (89) 


Is  found  from  (91)  and  is. 
1 


77  2 


3 


B 


*  r(§) 


9*9. 


m 


OB) 


(3/2\)~8  /3»(  0) 


o  s  9  <  em 


(100) 
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The  values  of  the  constants  are  completed  by  putting  in  the  range  < 


H(y) 


A'  =  -  |\  — exp l-A  ✓q  dy'Jdy 


from  (59).  On  collecting  our  results  together  we  have 
1  fa  H(y)  -Xj//q'dy' 


«M0)  =  r  I  —re 

J0  ** 


dy  ,  0_  <  6  <  - 


do)(0)  =  -  • 


Ai 


cco 


H(y)  { 

-  1?  Ki(T7)ds  ,  8 

a  a 


=  8„ 


and 


do>(0)  =  -- 


1  i 

3*  2®  H(y0) 


1  l  1 

^  d'cyo)2  q;4  ^U.iWj)  -  Jadi;) 


,  o  <6  <0, 


for  values  of  -q*  <  1  and  by  (105)  with  qw  replaced  by  q*  for  higher 

X  • 

We  can  approximate  to  (106)  by  putting 

1  hr  fy 

t)sKi(71)  ~  J—  e*’3  ,  where  17  =  A.)  \/q"  dy'  , 

3  V  2  Jo 


and  hence 


d«0>  -  -  — - - - -  [  ™  dy  .  «  = 

J.  ’• 


having  a  similar  form  to  that  given  in  the  range  9  <8  <  tt/2  . 


:  8  <  v/2 

(104) 

(105) 

(106) 

,  (107) 

values  of 

(106') 
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The  values  of  and  q'(y0)  can  be  obtained  from  (47)  and  If  we  ignore  deriva¬ 

tives  of  the  speed  of  sound,  noting  that  this  is  a  poor  approximation*  near  y  =  0  , 
we  find, 

q'  =  2\s|cok,|  (108) 


=  1  at  y  =  0 


=  I-3J  +  W 


where  Uj  and  5  are  evaluated  at  y  =  y0  . 

In  (109)  we  find,  on  approximating  to  q(y0)  =  0  ,  that 


q'(y0)  - 


and  from  (108) 


Hence,  with 


=  2X.9k|kJ 


(yO  ft  \  /& 

- •)  —  1  2ik,rdz, 

w  Vv 


(ill)  V 


as  given  by  (71),  we  can  find  expressions  for  Aw  ,  and  the  corresponding  values  of 
the  spectrum  functions  are 


_  4k?\2  dU.  3a e~'n  .  .  di),  pae~v‘  tt 

n(0;J5,,w)  =  — $(y,z;Jt/o)  — 1 - —  dy  — *■ - - —  dz  ,  9  <  9  <-  (112) 

dy  q*  dz  Q*  2 


n(.o;j<,cu) 


477kJ\2 


dUj  j>Ze~v  dU,  pa  e 

-  $(y,z;k/o)  — - t  -  dy  — - j - dz  , 

I  ~  dy  q*  dz  q* 


*  The  full  oxpres8ion  for  |qj  is  given  by 

| qw |  -  X2k 2 [1  -  <v2aV  -  (y-l)/2XV] 


but  note  our  previous  remarks  about  the  term  a"/a  • 
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We  find,  on  making  similar  approximations  to  those  used  above  in  evaluating  (81),  that 
the  integration  over  frequency  gives 


On  approximating  to  this  integral,  we  find  that 


h5(U)  1‘  L(Mro) 


(118) 


Just  as  in  the  previous  case  (83).  Thus,  although  we  have  taken  account  of  the  con¬ 
tribution  due  to  the  eddy  Mach  waves,  we  see  that,  according  to  the  approximations 
made,  they  do  not  contribute  more  than  those  eddies  which  travel  at  subsonic  speeds 
relative  to  the  wall. 


4,  CONCLUSIONS 

It  has  been  shown  that  theory  and  experiment  are  in  fairly  good  agreement  in  the 
prediction  of  wall  pressure  fluctuations  in  a  turbulent  Incompressible  boundary  layer. 

On  the  other  hand,  theory  and  experiment,  in  the  case  of  supersonic  flows,  show  some 
divergence,  although  both  appear  to  tend  to  a  similar  asymptotic  value  and  both 
demonstrate  the  presence  of  eddy  Mach  waves  above  a  certain  Mach  number.  However  the 
theory  does  not  show  any  marked  Increase  In  pressure  level  at  the  wall  due  to  eddy 
Mach  waves.  In  fact,  a  simple  extension  to  the  incompressible  flow  theory  is  shown  to 
give  similar  results  to  the  more  elaborate  theory  in  which  eddy  Mach  waves  are  approx¬ 
imately  taken  into  account.  This  result  is  perhaps  not  surprising  when  we  note  that 
the  dominant  region  associated  with  the  pressure  at  the  wall,  even  in  the  Incompressible 
case,  is  displaced  away  from  the  wall  and  the  level  of  pressure  is  roughly  constant 
over  an  appreciable  distance  normal  to  the  wall. 

It  is  clear  that  further  experimental  results  are  required  to  explain  thg 
differences  between  theory  and  experiment,  as  well  as  to  obtain  further  information 
in  the  case  of  flows  with  pressure  gradient  and  with  heat  transfer. 
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APPENDIX 


The  velocity  spectrum  function  for  u 

In  order  to  obtain  numerical  values  for  the  wall  pressure,  the  following  form  has 
been  chosen  for  the  velocity  spectrum  function  in  fixed  co-ordinates 

<*f>tf<xt  +  x')iJI2i32V 

- 1  “ 

2tt2  u0(i  +  I]**) 


$2j(x,,Xj  i  x';Js,,k}1£u)  = 


e'<lV+lX/» 


2 

0 


where  U0  is  the  eonvaction  speed  and  (sk!  ■  l*kj  +  i |k ®  ,  It  is  assumed  that  ls 
and  13  are  constant,  while  l2  -  x2  ,  thus  an  allowanoe  is  made  for  a  change  in 
scale  of  the  turbulence,  with  increase  in  distance  from  the  wall. 


The  Integration  over  all  values  of  co  leads  to 


$28(x2.x8  +  x»;k*8)  l\h2e’lV 

^  Juf(Xs  t  X')  77z(l  +  x|«J) 


and  its  Fourier  transform,  with  respect  to  *2  gives  the  two-dimensional  wave  number 
spectrum  function 


1k2x' 


$jj(x2,x2  +  x';k)  =  *22(x2,x2  +  x';^«2)e  z.  d*2 


_  e’** 


which  clearly  displays  qualitatively  the  correct  physical  properties  across  the 
boundary  layer,  even  though  it  fails  to  demonstrate  the  true  anisotropic  wave  number 
distribution. 


k /**.*./«  4- 4 


$22(X*,X2  +  x'ikj)  a 


uj(x2)  Ju |(X2  + 
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which  is  also  equal  to  $}4(x,,,  x?  +  x':  -  eu/U0)  in  agreement  with  Taylor’s  hypothesis. 
Finally  the  integration  over  all  kj  gives 


Rjj(0,  XjO;  0,  xs  +  x£,  0) 


ua(xa)ug(x„  +  x') 
ifuf(Xt)  )[uf(X4  +  x') 


a  torn  for  the  velnoit.y  nnrrelatlon  which  is  in  fair  agreement  with  experiment  over 
most  of  the  boundary  layer.  It  is  therefore  reasonable  to  expeot  that  our  calculated 
value  for  the  mean  pressure  which  is  based  on  this  assumed  form  for  $4J  win  be  in 
fair  agreement  with  its  exact  value  provided  very  high  wave  numbers  are  exoluded. 
Finally  it  is  worth  noting  that  we  have  a  complete  freedom  of  ohoice  with  respect  to 
the  longitudinal  and  transverse  soales,  and  l3  respectively,  as  well  as  the 
conveotlon  speed  U0  ,  except  that  in  the  latter  case  we  have  assumed  that  it  is 
independent  both  of  the  distance  from  the  wall  and  of  the  wave  number  in  the  turbulenoe. 
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DISCUSSION 


G.M.  Corcos 

Professor  Lilley  has  based  much  of  his  talk  upon  the  assumption  that  the  pressure 
field  was  caused  primarily  by  the  quasi  linear  momentum  flux  which  results  from  the 
interaction  of  the  mean  stream  rate  and  one  turbulent  strain  rate.  He  and  I  have  been 
having  a  disagreement  on  this  point  which  could  hardly  be  more  complete  or  more 
friendly.  He  has  advanoed  many  arguments  in  favor  of  his  hypothesis  originally  pro¬ 
posed  by  Vralchan  and  this  Is  not  the  place  to  review  eaoh  one  in  detail  as  I  had 
originally  planned  to  (slide  down).  Fortunately  the  discussion  has  shifted  from  the 
somewhat  speculative  plane  of  the  appropriateness  of  certain  modes  of  things  we  don’ t 
know  to  that  of  hard  faots.  By  this  I  mean  that  the  excellent  measurements  of 
Wooldridge  and  Wlllmarth  allow  the  direct  evaluation  of  the  contribution  of  this  one 
inertia  term  to  the  pressure  field  without  recourse  to  models  or  assumptions.  One 
then  can  compare  these  contributions  to  the  characteristics  of  the  observed  pressure 
field.  Now  Professor  Lilley  has  almost  done  that  by  trying  to  fit  a  pressure-normal 
velocity  correlation  model  to  these  experimental  data.  This  later  approach  I  have  not 
had  time  to  check  in  detail,  since  I  only  got  the  report  yesterday.  But  there  seem 
to  be  Important  discrepancies  between  the  model  and  the  measurement  near  the  wall. 

If  there  were  not,  which  is  still  possible,  Ills  computation  should  agree  with  a 
detailed  numerical  computation  I  have  performed  without  a  model  and  therefore  without 
assumptions,  unless  I  made  many  numerical  mistakes.  Briefly  the  results  of  these  two 
computations  are  different  as  day  and  night.  Professor  Lilley  finds  that  the  M-T  or 
quasi  linear  term  dominates,  X  find  it  does  not  dominate.  Professor  Lilley  finds 
that  this  term  makes  its  contribution  around  y*  =  2000  or  y/8*  =  1.6,  I  find  it 
makes  its  peak  contribution  around  y  =  y*  =  350  -  400  ,  y/S*  .2  .  In  addition  I 
find  that  the  auto-correlation  of  the  pressure  contributed  by  that  spectrum  is  very 
different  from  the  observed  correlation  -  contributing  essentially  only  the  very  high 
frequencies  of  the  spectrum.  The  issue  can  clearly  and  relatively  easily  be  resolved. 
We  shall  first  compare  the  model  to  be  usBd  with  the  data  in  detail  and  if  they  do 
agree,  then  we  shall  look  for  numerical  mistakes,  I  intend  to  describe  briefly  this 
afternoon  how  the  numerical  computation  of  that  contribution  is  made  "811x1  to  give  a  few 
details  of  its  results.  I  would  like  to  note  that  even  if  allowance. is  made  for  this 
rather  Important  assumption  with  which  I  disagree,  Professor  Lilley’ s  talk  was 
thoroughly  enjoyable  and  his  paperwork  on  compressible  flow  is  quite  interesting. 

Author's  reply 

I  would  disagree  that  I  had  based  my  work  primarily  on  the  importance  of  the  mean 
shear-turbulence  interaction.  ' 


J.E.  Ffouics  Williams 

I  was  particularly  interested  to  hear  of  Professor  Lilley’ s  work  on  the  wall 
pressure  fluctuations  under  compressible  flow  conditions.  I  would  like  to  ask 
Professor  Lilley  if  he  has  given  any  thought  to  one  aspect  of  his  theory  about  which 
I  feel  uneasy.  That  is  concerned  with  the  problem  of  whether  or  not  his  assumed  model 
presents  a  well  posed  mathematical  problem.  The  model  is  that  of  an  infinite  plane 
supporting  a  homogeneous  boundary  layer.  The  turbulence  will  generate  sound 


and  we  are  interested  in  that  sound.  The  situation  is  such  that  it  gives  rise  to  an 
analogue  Of  Okert' s  paradox  for  the  pressure  at  any  one  point  is  infinite.  The 
infinity  is  avoided  in  Professor  Lilley’s  equations  by  restricting  the  analysis  to 
frequencies  given  by  KjU  .  But  the  total  pressure  is  nevertheless  infinite  because 
finite  contributions  come  from  all  parts  on  the  plane,  i.e.  f dS/r*  is  singular.  Is 
the  analysis  still  valid  even  though  it  has  to  cope  with  this  infinite  contribution? 

,4u tho r ' s  reply 

I  have  looked  at  this  problem  and  I  felt  it  was  sufficient  to  assume  that  the  plate 
was  not  necessarily  Infinite  but  large.  It  would  appear  that  the  results  I  have 
obtained  for  the  pressure  covariance  do  in  fact  converge. 
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